Scale dependence of branching in arterial and bronchial trees 
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Although models of branching in arterial and bronchial trees often predict a dependence of bi- 
furcation parameters on the scale of the bifurcating vessels, direct verifications of this dependence 
by comparison with data are uncommon. We compare measurements of bifurcation parameters of 
airways and arterial trees of different mammals as a function of scale to general features predicted 
by theoretical models based on minimization of pumping power and network volume. We find that 
the size dependence is more complex than existing theories based solely on energy and volume 
minimization, and suggest additional factors that may govern the branching at different scales. 
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The potential factors determining the parameters at 
bifurcations in the arterial and bronchial trees of mam- 
mals have lately received attention. Most notably, these 
local bifurcation characteristics have been proposed to 
determine allometric scaling laws for biologically impor- 
tant variables (see 0,0 and 0). It has been suggested 
that arterial and bronchial trees follow a pattern that 
minimizes a combination of pumping power and the to- 
tal volume of the tree (see for example If the 
flow is assumed to be Poiseuille (resistance proportional 
to the inverse fourth power of the vessel radius) on all 
scales, self-similar branching results, i.e., the branching 
parameters do not depend on the scale. However, it has 
been argued Qltliat, due to the fact that the blood flow 
is pulsatile Oj , arterial trees determined by optimiza- 
tion have branching parameters that depend on the scale 
(the form of this dependence is discussed below). As 
a direct test of these theories and their assumptions, we 
compare in this Letter experimental measurements of the 
scale dependence of the branching parameters with the 
predictions of theoretical models. Many studies of bifur- 
cation parameters rely on ordering schemes [Tol | which 
do not directly consider the dependence of the bifurca- 
tion parameters on the scale, and thus are not appropri- 
ate for our particular purpose. We will explicitly study 
the dependence of the radii of the bifurcated vessels on 
the radius of the parent vessel, thus allowing a more di- 
rect comparison with theoretical models. We show that 
important aspects of previous observations |ll| for the 
human bronchial tree, not predicted by the theoretical 
models, are more general than has been previously no- 
ticed. In particular, we find qualitatively similar behavior 
for dog arterial and bronchial trees, and lamb fetal and 
neonatal arterial trees. Our results suggest that, in addi- 
tion to minimizing pumping power and network volume, 
there might be other important factors determining the 
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branching in the arterial and bronchial trees. 

Arterial and bronchial trees start with a tree root (a 
parent vessel) that divides into two (daughter) vessels, 
and these in turn subdivide in a similar fashion, contin- 
uing until a certain approximate value of the radius is 
attained (in the case of arterial trees, the capillary ra- 
dius). At a particular bifurcation, wc let ro denote the 
radius of the parent vessel and r\ , denote the radii of 
the daughter vessels. Variables often used to character- 
ize the bifurcations are the area ratio and the symme- 
try index, defined respectively by a = (r 2 + r|)/?"Q and 
a — T2/ri, where r2 < ri [T^. These two variables com- 
pletely determine r\ and r-i in terms of r$. The branching 
exponent x defined for ro > r\ > r-2 by rfi = rf +rf is also 
often used in modeling and theoretical discussions |l3j- 
[T^ | . Most predictions of theoretical models are stated in 
terms of this exponent. However, we find this exponent 
inconvenient for our purposes since it diverges as r\ or 
r-2, approach ro, as occurs frequently with experimental 
data. We will therefore use the area ratio in our analysis. 
The value of the exponent x can be linked to the change 
of certain quantities at bifurcations. An exponent x = 2 
indicates conservation of total cross sectional area (a = 1) 
and a value x = 3 corresponds to conservation of laminar 
shear stress at a bifurcation. The area ratio a is related to 
x and the symmetry index a by a = (1 + a 2 )(l +a x )~~ . 
This is a monotonic function of x which does not change 
much as the value of a is varied between 0.5 and 1. We 
will therefore regard the general features of the area ra- 
tio as a function of the scale as a counterpart of a similar 
qualitative behavior of the exponent x. 

Murray 3 proposed that the arterial tree is con- 
structed in such a way as to minimize a linear combi- 
nation of the pumping power and the volume of the net- 
work. This leads, if Poiseuille flow is assumed through- 
out the tree, to Murray's Law, which states that x = 3. 
Murray's Law has been useful as a reference point regard- 
ing the small blood vessels and is sometimes mentioned 
in discussions regarding the bronchial tree as well pflj . 
However, deviations have been apparent for some time 
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(see 0|" an d references therein). The exponent x 
has been observed to be generally smaller than the pro- 
posed value of 3. 

Recently, West et al. [l[ included the effect of the 
blood flow being pulsatile, and found that the minimiza- 
tion scheme carried on by Murray leads in this case to 
a transition scenario, with x close to 2 for the larger 
vessels and to 3 for the smaller vessels (see, however, 
Ref. 0]). This modification predicts a smooth transition 
from the larger to the smaller vessels, as suggested also by 
Uylings 5] on other grounds. Other modifications have 
been proposed predicting an exponent between 2 and 3 
for the smaller vessels 00- I n order to test these the- 
ories we study data for the bifurcation parameter a as a 
function of the parent vessel radius. Zamir [l2T | presents 
data for the branching parameters as a function of the 
radius for the right coronary artery of a human heart. 
The large scatter in the data, however, makes it diffi- 
cult to observe the dependence of the mean values on the 
scale. Phillips and Kaye studied the dependence of 
averaged branching parameters on the radius of the ves- 
sels in the bronchial tree of four mammalian species. We 
show that some important aspects of their observations 



for the human bronchial tree (to be discussed) are more 
general than has been previously noticed. 

Data sets. We study morphological data from the 
bronchial tree of a human and 2 dogs obtained by Raabe 
et al. [24[, a human bronchial tree from Horsfield et al. 
plj . dog lung arterial trees from Dawson et al. 0] and 
neonatal and fetal lamb lung arterial trees from Bennett 
et al. 0] (see Table 1). For each data set, we compute 
the average value for the area ratio (a) as a function of 
the parent vessel radius r as follows. For a given parent 
vessel radius r we construct a sample of all bifurcations 
(ro,ri,ra) which satisfy rlO~ c < ro < rlO e , thus making 
uniform intervals in logr, and take (a) evaluated at r to 
be the average of a over this finite sample. We choose 
the logarithmic interval width 2e to be 0.05. This value 
was chosen so that the statistical noise is filtered but 
the large scale features of (a) as a function of r are not 
smoothed out. For the largest vessels and sometimes for 
the very small vessels, the number of bifurcations mea- 
sured is sparse, especially when only one lung was mea- 
sured. We reject the values obtained when the sample 
for a given range of radii contains less than 4 members. 



Data set 


Mammal 


Tree 


Body mass (kg) 


Radius range (cm) 


Bifurcations 


Horsfield et al. 


Human 1 


Bronchial 




0.07 - 1.6 


3088 


Raabe et al. 


Human 2 


Bronchial 


81 


0.04 - 2.01 


4336 


Raabe et al. 


Dog 1 


Bronchial 


11.6 


0.02 - 1.8 


8032 


Raabe et al. 


Dog 2 


Bronchial 


10.3 


0.03 - 1.61 


2057 


Dawson et al. 


Dogs (32) 


Lung Arterial 


19.7 ± 2.1 


0.003 - 0.76 


919 


Bennett et al. 


Lamb (fetal) 


Lung Arterial 




0.001 - 0.70 


10970 


Bennett et al. 


Lamb (neonatal) 


Lung Arterial 




0.011 - 1.23 


846 



TABLE I: Characteristics of the different data sets used. A 
"-" in the body mass entry indicates that this information 
was not available 

The results of this procedure are shown in Figure ^ 
In order to allow comparisons with the values of the ex- 
ponent x predicted by theoretical models, we display the 
values of a corresponding to thick line on the 

right axis. The corresponding value of the area ratio a 
for x = 3 depends on the symmetry index a, and thus we 
marked the range of values of a corresponding to x — 3 
for values of a ranging from 0.5 to 1. The error bars 
correspond to one standard deviation as computed from 
the propagated measurement error taken from the preci- 
sion to which the data were reported, and the statistical 
uncertainty due to a finite sample. 

In Figs. I2a) and (b) we show (a) as a function of the 
parent radius r for the human and dog bronchial trees 
respectively, obtained as described above from the Raabe 
et al. and Horsfield et al. data sets. In Fig. QIc) we 
display the area ratio from the dog pulmonary arterial 
tree, and fetal and neonatal lamb arterial trees from the 



data of Dawson et al. and Bennett et al., respectively. 
The behavior for the largest vessels is not shown in the 
plots because of lack of good statistics. It has been shown 
elsewhere |53 that the bifurcations of the largest vessels 
in the arterial tree generally have an area ratio close to 
1. 

We first discuss the results for the dog and human 
bronchial trees [Figs-^a) and (b)]. The behavior of (a) 
as a function of r is not a simple one: it is not constant, 
as it would be if the branching were self-similar, nor does 
it have a smooth transition from a lower value to a higher 
value as r decreases. Remarkably, the same qualitative 
behavior is observed in both graphs. This behavior con- 
sists of the following pattern: for the very small and close 
to terminal vessels, the area ratio is high. As r increases, 
(a) decreases, reaching a local minimum at r ~ 1 — 2 mm. 
Then the area ratio grows, and, for Human 1 (boxes in 
[Fig. Ola)]) and the dogs [Fig. ^b)], one can observe it 
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(a) Human airways 
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(b) Dog airways 
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■ (c) Lung arteries 
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FIG. 1: Average area ratio (a) as a function of the parent 
airway radius r for (a) Horsfield et al. Human 1 (boxes) and 
Raabe et al. Human 2 (diamonds), (b) Raabe et al. Dog 1 
(triangles) and Dog 2 (stars) (see Table 1), and (c) Dawson 
et al. dogs lung arterial tree (triangles) and Bennett et al. 
neonatal (stars) and fetal (diamonds) lambs lung arterial tree. 



finally decreasing again. In quantitative terms, we make 
the following observations. For the humans, among the 
largest vessels, the value of the area ratio (a) peaks higher 
than that corresponding to a value of x = 3, even when 
the error bars are taken into account. For smaller ves- 
sels, (a) decreases as mentioned before, attaining values 
below those equivalent to i = 3. For the smallest ves- 
sels, the area ratio takes values substantially above those 
corresponding to x = 3. 



We observe a qualitatively similar situation for the 
lung arteries [Fig. [He)]. For each data set, there is a 
local minimum of the area ratio at moderately small ves- 
sel radius, and the area ratio increases sharply as the ra- 
dius decreases further. The fetal and neonatal lamb arte- 
rial trees (diamonds and stars, respectively) have smaller 
overall values of the area ratio, as noted before In 
the dog arterial tree (triangles) the different regimes are 
particularly marked: area ratio close to 1 for the largest 
vessels, a peak for the moderately large vessels, a smaller 
value of the area ratio as the radius decreases further, and 
larger values of the area ratio near the terminal vessels. 
The presence of this peak for moderately large vessels 
was noted by Phillips et al. ^3] for the bronchial tree 
of the human in Raabe et al. data set; we see that this 
behavior may be more general. 

We also studied data for the airways of two rats and a 
hamster from Raabe et al. [IH , but the range of vessel 
radii was smaller and the results were not conclusive. 
Data from the heart arteries of a pig from Kassab et 
al. |2fj also did not show a clear dependence of area 
ratio on vessel radius. While not inconsistent with our 
observations above, they were not inconsistent with the 
hypothesis of a constant area ratio. 

The above discussion is qualitative, and it is natural 
to ask whether it can be supported by a rigorous statis- 
tical test. We believe, however, that due to the limited 
amount of samples available and to the inherent differ- 
ences between each of the sampled individuals, we do not 
have sufficient basis for formulating an alternative model 
to those we have discussed. We hope that the graphical 
information we have presented will encourage experimen- 
talists to measure and make available more data. 

In summary, a simple transition of the area ratio from 
values close to 1 for the large vessels to higher values for 
the smaller vessels, as predicted by some optimization 
models, is not supported by the data we studied. The 
assumption that an optimization principle is involved in 
the design of the arterial and bronchial trees is natural, 
as evolution would tend to eliminate designs that are not 
adapted to meet the demands of the organisms. We be- 
lieve, however, that these demands are not restricted to 
energy efficiency (as assumed in most theories) and that 
different factors compete, resulting in the complex be- 
havior of the area ratio. As has been remarked recently 
|28| , a bronchial tree designed solely with the purpose of 
minimizing pumping power and network volume could be 
dangerously sensitive to deviations from the optimal bi- 
furcation values, and higher area ratios could be a safety 
margin against the effect of these variations. Other fac- 
tors possibly contributing to make the area ratio smaller 
for the largest vessels are impedance matching and mini- 
mization of the resistance for pulsatile flow [l] . A higher- 
value of the area ratio as the radius decreases could in- 
dicate a tendency to relieve the high pressure and shear 
stress present in the large vessels. The large value of the 
area ratio near the terminal vessels could allow the fluid 
to be slowed down as it reaches the capillaries or alveoli 
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Higher area ratio 



Minimize steady flow resistance 



Relieve high 

Slow fluid down pressure / shear stress 



Terminal Large 
vessels vessels 



Minimize 

pulsatile resistance 

Match 

impedances 



Minimize volume / average time 

Lower area ratio 

FIG. 2: Schematic representation of possible factors affecting 
the area ratio on different scales. The factors above the thick 
horizontal line could favor higher area ratios, and those be- 
low the horizontal line favor a lower area ratio (our vertical 
order is otherwise arbitrary). A compromise between these 
competing factors could produce transitions as seen in Fig. 1. 



(e.g., to allow fuller transfer of the transported nutrient 



to the body). On the other hand, minimizing the volume 
or the average circulation time favors a smaller area ra- 
tio on all scales. Also, the effects of finite vessel length 
and vessel curvature could potentially affect the values 
of the optimal branching parameters. This discussion is 
summarized in Fig. EI We emphasize that these are only 
possible factors, and that this list is not meant to be ex- 
haustive. In conclusion, we have presented evidence of 
a consistent and so far unexplained behavior in diverse 
arterial and bronchial trees. To explain these observa- 
tions it might be necessary to couple existing optimality 
principles to other requirements. The precise formula- 
tion and understanding of such constraints is a nontrivial 
and important problem. The consequences for allomctric 
scaling are also important, as the branching pattern of 
the arterial tree has been argued |l| to be determinant 
of the scaling of many variables. We believe it will be 
fruitful to carry on a similar analysis on other data. 
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